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Abstract. We explore the feasibility of using the Moon as a detector of extremely 
high energy (> 10^^ eV) cosmic rays and neutrinos. The idea is to use the existing 
radiotelescopes on Earth to look for short pulses of Cherenkov radiation in the GHz 
range emitted by showers induced just below the surface of the Moon when cosmic 
rays or neutrinos strike it. We estimate the energy threshold of the technique and the 
effective aperture and volume of the Moon for this detection. We apply our calculation 
to obtain the expected event rates from the observed cosmic ray flux and several 
representative theoretical neutrino fluxes. 



INTRODUCTION 

The observation of atmospheric showers with total energy above the so-called 
GZK cutoff at energies above ~ 5 x 10^^ eV, is one of the most puzzling mysteries 
of the emerging field of astroparticle physics [1] . The nature of the primaries is still 
unknown as well as their origin and the mechanisms involved in their acceleration 
to such energies. The puzzle is even more intriguing because if the primaries are 
protons, photons or nuclei they should be produced within distances of a few tens 
of megaparsecs from us, otherwise their interactions with the 2.7 K photons consti- 
tuting the Cosmic Microwave Background (CMB) as well as with the infrared and 
radiobackgrounds, should greatly reduce their energy. Since no effective sources 
capable of accelerating them have been identified within that distance, their very 
presence is contrary to the expectations. Neutrinos may play a fundamental role in 
solving this puzzle. They may be produced in interactions of the primary cosmic 
rays with the matter and radiation in the source or in the matter and radiation 
they find along their paths to Earth. The so-called GZK neutrinos produced in 
interactions of cosmic rays with the CMB are almost guaranteed, since both the 
projectile and the target are known to exist [2,3]. Sources in which ultra high en- 
ergy neutrinos may be produced include AGN's and GRB's [4] where accelerated 



particles interact with matter or radiation. Alternative sources are the annihila- 
tion of topological defects which can be produced in phase transitions in the early 
universe where extremely high densities of energy are trapped [5]. 

Some present and future detectors of ultra high energy cosmic rays [1,6-8] and 
neutrinos [4,9,10] will soon remedy the lack of statistics [11]. The small observed 
(predicted) fluxes of ultra high energy cosmic rays (neutrinos) call for huge effec- 
tive volumes (measured in km^). Several detection methods and techniques are 
being employed and/or studied to achieve the required volumes, some of them have 
been discussed in this workshop [12]. Among them the radio technique is one of 
the most promising alternatives for neutrino and possibly cosmic ray detection at 
ultra high energies (10^^ eV and above). The technique aims at detecting coherent 
Cherenkov radiation in the MHz-GHz range from the excess of electrons in showers 
initiated by photons or electrons. Several detectors are being planned [12] or they 
are in the early stages of construction [13]. Due to its excellent radio frequency 
wave propagation properties, Antarctic ice is being considered as a medium where 
antennas are being placed to monitor the potential radio signals. Other media such 
as salt [14] and sand are also being studied. As an interesting alternative, in 1989 
Zheleznykh et al. [15] proposed to detect showers initiated by cosmic ray and neu- 
trinos by measuring coherent Cherenkov radiation emitted just below the surface 
of the Moon when high energy cosmic rays or neutrinos strike it. Two groups have 
looked for those signals using existing radiotelescopes on Earth with no positive 
detection during the time they have pointed the instruments to the Moon [16,17]. 

Here we investigate the feasibility of detecting extremely high energy particles 
interacting in the Moon. We discuss the characteristics of the radio signals and 
we stress how they influence the aperture of the Moon as a high energy particle 
detector. We calculate the energy threshold for the technique as a function of the 
detector parameters and the signal's main features. Using several representative 
neutrino fluxes, we estimate the expected event rates and we briefly discuss where, 
on the surface of the Moon, should the radiotelescopes be pointed at to be able to 
detect them. We also discuss cosmic ray detection. This work updates previous 
estimates in which some approximations to extremely high energy shower develop- 
ment in the lunar rock as well as to the associated radiosignals were used [18]. We 
also take into account the transmissivity of the signals through the Moon's surface 
and some other geometrical issues which were not fully considered in our previous 
estimates. 

RADIOPULSES IN THE MOON 

Neutrino interactions produce different types of showers depending on the neu- 
trino flavor and on whether the interaction is mediated by a W^ or a Z boson. In 
deep inelastic scattering (DIS) charged current (CC) interactions of electron neu- 
trinos (z/g), the electron produced in the lepton vertex initiates an electromagnetic 
shower of energy (1 — y)Ej,, where Ej, is the neutrino energy in the laboratory frame 



and y is the fraction of energy transferred to the struck nucleon in the hadronic 
vertex. The debris of the nucleon initiate a hadronic shower which is superimposed 

to the electromagnetic one producing a "mixed shower" . Hadronic showers are also 
initiated in both charged and neutral current (NC) DIS neutrino interactions by 
the hadrons created in the fragmentation of the nuclear debris. At shower energies 
above -Blpm ~ 4 x 10^^ eV (400 TeV) in the lunar regohth, electromagnetic showers 
are affected by the LPM effect [19]. When the energy of a photon or electron in the 
shower is above -Elpm the interaction distance becomes comparable to the inter- 
atomic spacing and collective atomic and molecular effects affect the static electric 
field responsible for the interaction. The result is a reduction in both total cross 
sections with lab energy [E) which drop like E~^-^ above ~ Elpm- As a result 
PeV - EeV photon or electron induced showers are considerably larger than at TeV 
energies [20]. In hadronic showers this effect is mitigated because the main source 
of photons is the decay of 7r°'s which at energies above E^^o = 7 x 10^^ eV (7 PeV) 
is suppressed due to the dominance of 7r° interaction [21]. 

A helpful way of getting insight into the pattern of the electric field emitted by 
the excess of electrons in the shower, is to think of a shower as a slit illuminated by 
radio- waves [22]. The diffraction pattern can be obtained as the Fourier transform 
of the slit, however due to the nature of the Cherenkov radiation the emission 
has a central peak with its maximum at the Cherenkov angle instead of in the 
direction perpendicular to the slit. The radiation is coherent when the wavelength 
is larger than the physical dimensions of the shower. In this situation the radiated 
electric field becomes proportional to the excess charge and the power in radiowaves 
scales with the square of the shower energy. This effect was predicted in the 1960's 
by Askary'an [23] and has been observed recently at SLAC in experiments where 
an intense beam of GeV photons generates a shower inside a sand target [24]. 
Given this interpretation it is easy to understand that the LPM effect is going to 
reduce the angular width of the Cherenkov peak (A6') in electromagnetic showers 
as shower energy increases. Following the same reasoning, the diffraction pattern 
from hadronic showers is expected to be less affected by the LPM, exhibiting a mild 
decrease of A6' with energy [21]. 

The spectrum of the electric field also increases linearly with frequency up to 
a maximum frequency which is essentially determined by the lateral structure of 
the shower at the Cherenkov angle, the wider the shower the smaller the corre- 
sponding maximum frequency. At angles away from the Cherenkov angle, the 
maximum frequency depends only on the longitudinal dimension of the shower. 
We have performed simulations of electromagnetic showers in the lunar regolith as- 
sumed homogeneous (density p ~ 1.6 g cm~^). We have adapted the ZHS code [25] 
originally conceived for simulations in ice, changing the relevant parameters such 
as density, radiation length, atomic number and refraction index. The frequency 
spectrum as well as the angular distribution of the pulses emitted by a 100 TeV 
electromagnetic shower in the Moon are shown in Fig.l. The results of the simu- 
lations are in very good agreement with a simple scaling of the simulations from 
ice to the Moon [18]. With this in mind we have obtained the electric field from 



electromagnetic showers at the highest energies, by scahng the spectral features in 
ice. The same scahng has been apphed to estimate the radiopulses from hadronic 
showers. The absolute value of the electric field has been normalized according 
to the amount of electromagnetic energy present in the hadronic shower, and the 
width of the Cherenkov peak, which is inversely proportional to the longitudinal 
dimension of the shower, is scaled with the radiation length and the density of the 
medium [18]. Simulation work is in progress in this direction and will be published 
elsewhere. For more details see references [20,21]. 
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FIGURE 1. Left panel: Results of simulations of the angular behavior of the electric field 
emitted by a 100 TeV electromagnetic shower in the lunar regolith. 9 is the observation angle 
with respect to the shower axis. Two frequencies are shown. Right panel: Frequency spectrum 
of the electric field at two observation angles. 




APERTURE OF THE MOON FOR NEUTRINO AND 
COSMIC RAY DETECTION: ENERGY THRESHOLD 

One important parameter for assessing the possibilities of using the Moon as 
a high energy particle detector is the energy threshold for shower observation by 
radiotelescopes on Earth. This can be estimated comparing the radio signal emitted 
by an electromagnetic shower at the Cherenkov peak and the noise expected by 
other processes in the radiotelescope. Here we only take into account thermal noise 
for which we use two representative values, namely la and Gcr of the flux density 
given by the standard expression for a radio antenna: 



Here /cb is Boltzmann's constant, Tgys is the noise temperature of the radio detection 
system, At and Au are respectively the duration of the pulse and the bandwidth of 

the detection system around a central frequency uq- A^g is the effective area of the 
antenna. We have estimated the energy threshold for the NASA/JPL Goldstone 
Deep Space Station 14 (70 m diameter dish) radio telescope [26] with which the 
measurements in [17] were performed. The la thermal noise level for this system 
is ~ 400 Jy ^ and the Qa is 2,400 Jy. The flux density on Earth emitted by an 
electromagnetic shower at the Cherenkov peak is obtained integrating the frequency 
spectrum of the electric field RE(i/, R, 6'cher), where R is the distance to the shower, 
around uc- 
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REarth-»Moon is the distaucc from Earth to the Moon. The result of this calculation 
for electromagnetic showers is shown in Fig. 2. Reading from the plot one can 
estimate the energy threshold of a electromagnetic shower to be between Et^ ~ 
10^^ — 10^*^ eV depending on the noise level, the central frequency and assuming a 
bandwidth of ~ 0.1 i/c- To translate the electromagnetic shower energy threshold 
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FIGURE 2. Left panel: Flux density at Earth emitted by an electromagnetic shower of energy 
-E'eiectrom. for two Central frequencies and a bandwidth Au = O.lvc- Also shown is the flux density 
of the thermal noise in NASA/JPL Goldstone DSS 14 radiotelescope. Right panel: Transmissivity 
(Power transmitted/Powerincident) of radio wavcs at the Moon-vacuum interface as a function of the 
incidence angle with respect to the perpendicular to the local surface of the Moon. 



into a neutrino or cosmic ray energy threshold, one has to keep in mind that at 



1) IJy = 10-26 W m-2 Hz-i 



extremely high energy the electromagnetic energy content of a hadronic shower in 
the lunar regolith is ~ 90% of the total shower energy, slowly increasing with it [21] 

and the typical value of y is ~ 0.2 [27]. This implies that the energy threshold for 
detecting a NC neutrino DIS interaction is roughly E^y ~ 5£^th whereas it is ~ £^th 
for detection of CR or for a CC i^e interaction. 

Aperture 

We have studied in detail the geometry of the problem in order to calculate the 
aperture (effective area x solid angle) of the Moon for neutrino as well as cosmic 
ray detection. 

For a cosmic ray or a neutrino to be detected it has to interact producing a 
shower within a distance to the surface of the Moon equal to the typical absorption 
length of radiowaves inside the lunar regolith: X^^^° ~ 15 m(l GHz/z/). The emitted 
radiowave then reaches the surface of the Moon without significant attenuation and 
has to point to the Earth after refraction in the Moon- vacuum interface. We restrict 
ourselves to radio emission from the directions between the Cherenkov angle and 
6c ^^6, where most of the power is concentrated. These two considerations allowed 
us to obtain the potentially detectable incident neutrino and cosmic ray directions. 
Then we disregarded those directions for which either the cosmic rays or neutrinos 
are absorbed inside the Moon before reaching its surface^. It requires a httle bit 
of thinking and the aid of Fig. 3 to realize that this limits the potential interaction 
points on the surface of the Moon to a fairly narrow rim for cosmic rays and to a 
wider rim for neutrinos. It also restricts the observed cosmic rays to those hitting 
the Moon with their tracks almost parallel to its surface. We finally eliminated those 
directions (and energies) for which the signal, although being "geometrically able" 
to reach Earth, is below the noise flux density. In this respect the transmissivity 
properties of the Moon-vaccumm interface play a fundamental role. For a shower 
propagating inside the Moon parallel to the local surface, the radiation emitted 
at the Cherenkov angle suffers total internal reflection since the Cherenkov angle 
is the complementary of the total internal reflection angle. However detection is 
still possible because, as we mentioned above, the radiation is emitted in a cone 
of width A9 and the transmissivity of the interface increases to almost 1 just a 
few degrees off the total internal reflection angle as can be seen in Fig. 2. Our 
results on the aperture are shown in Fig.4. It is interesting to notice that for 
hadronic showers once the energy threshold is reached, the aperture varies slowly 
with shower energy, the reason for this being that the effective solid angle is not 
restricted by an energy decrease of caused by the LPM effect. This is not the 
case for electromagnetic showers where the decrease of the aperture with energy 
is induced by a corresponding reduction in due to the LPM. The inner panel 
in Fig.4 represents the visible face of the Moon and gives an idea of the size of 



^) The interaction length Lint"^ of a neutrino is smaller than the diameter of the Moon when 
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FIGURE 3. Main elements entering in the geometry of radio detection in the Moon. All the 
cosmic ray and neutrino directions entering in the depicted effective solid angle are "geometrically 
allowed" , although many of them are disregarded due to absorption inside the Moon or threshold 
effects (see text). 

the rim wliere cosmic ray and neutrino detection might be expected at 10^^ eV. 
When combining the aperture with A^^g'" acceptances in excess of 1000 km^ sr at 
El, ~ 10^° eV might be achieved. 



EVENT RATE ESTIMATES 

Cosmic ray event rates are calculated in a straightforward manner by convoluting 

the experimentally observed flux and the aperture of the Moon. For neutrinos the 
convolution involves the theoretically expected neutrino flux, the neutrino cross 
section and the effective volume of the detector: 

A^CR=/ ^^S^^'dEcK (yr-i); for cosmic rays, (3) 

r d^ \radio 

7V.= / :7^S^^-^dE, (yr-^); for neutrinos, (4) 

where S^^' ^ is the aperture of the detector. 

In Fig. 5 we show the event rates for the conservative case where the noise level is 
equal to Gcr of the thermal noise. We have used some neutrino fluxes representative 

of the different current models and theoretical ideas. For the cosmic ray event 
rate, we have used a simple parameterization of the experimentally observed flux 
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FIGURE 4. Aperture of the Moon for cosmic ray and neutrino detection. For neutrino initiated 
mixed showers the aperture is the sum of the sohd and dashed curves (not drawn). The inset 
represents the surface of the Moon. The hues correspond to the same hnes in the legend. The 
region between each line and the outermost solid line is the surface where cosmic ray or neutrino 
detection is expected at = 10^^ eV. 



which we have extended, perhaps conservatively, to the highest energy observed 
event (3 x 10^'' eV) using a differential slope 7 ~ —2.7. It is clear from the figure 
that cosmic ray events are going to dominate the event rate near the rim of the 
Moon outnumbering the neutrino events at least for the representative neutrino 
fluxes we have chosen. It is then advisable to point the radiotelescopes towards the 
center of the Moon to be able to observe neutrino events. Due to their high angular 
resolution, the radiotelescopes such as Goldstone DSS 14 or Parkes in Australia, 
are able to observe roughly 10% of the surface of the Moon at once, hence one 
should multiply the numbers in Fig. 5 by a factor of 0.1, to get an estimate of the 
sensitivity of the existing instruments and the observation time required to collect 
a certain number of events. Our estimates indicate that Goldstone DSS 14 should 
roughly expect 1 cosmic ray event every 80 hours of observation. 



IMPROVEMENTS AND CONCLUSIONS 

There are still a few issues that should be explored. The roughness of the surface 
of the Moon should be included in the calculation. Its wrinkles may facilitate 
the detection of signals emitted at the Cherenkov angle since, in this case, the 
particle tracks skimming the regolith surface are in general not parallel to the 
local surface which creates a problem of total internal reflection. Since showers 
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FIGURE 5. Cosmic ray and i^e CC + Ve NC + CC+NC events in the Moon expected in the 
NASA/JPL Goldstone DSS 14 radiotelescope. The events are shown as a function of the cosine 
of the latitude of the points on the surface of the Moon, namely cos(6'iat) = 1 corresponds to the 
center of the Moon and cos(^iat) = to the outermost rim. One should understand this plot as 
having azimuthal symmetry around cos(0iat) = 1- The numbers accompanying the legends are 
the total events per year. The numbers in parenthesis represent the neutrino events per year from 
the latitudes where cosmic ray events are not expected. For cosmic rays the number of events 
below 3 X 10^° eV is shown. 



should be produced near the surface of the Moon, one should worry about near- 
field effects in the Cherenkov emission [30] . Their main consequences are a decrease 
of the normalization of the electric field at the Cherenkov peak accompanied by 
an increase of the width of the Cherenkov cone [31]. These two effects act in 
opposite directions with respect to the behavior of the event rate. A decrease in 
will increase the effective sohd angle while the decrease in the normalization 
increases the energy threshold. There are also some issues dealing with showers 
produced at a distance to the surface of the Moon smaller than a radiation length 
so that Cherenkov radiation doesn't have enough distance to form. Variations of the 
index of refraction near the surface of the Moon have been measured and should also 
be included in the calculation of the geometry. Work is in progress to estimate the 
importance of all these potential issues and the results will be published elsewhere. 

Although some room for improvement docs exist, we have identified and correctly 
taken into account many of the elements entering in this difficult problem, updating 
previous estimates [18]. The actual potential of the Moon as a detector of high 
energy cosmic rays and neutrinos depends on a comprehensive study of all the 
elements indicated in this paper and possibly others. In any case, our results are 
encouraging and should trigger more theoretical work on this exciting possibility. 
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